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A Langevin approach to understand the noise of microwave devices is presented. The 
device is represented by its equivalent circuit with the internal noise sources included as 
stochastic processes. From the circuit network analysis, a stochastic integral equation 
for the output voltage is derived and from its power spectrum the noise figure as a func- 
tion of the operating frequency is obtained. The theoretical results have been compared 
with experimental data obtained by the characterization of an HEMT transistor series 
(NE20283A, by NEC) from 6 to 18 GHz at a low noise bias point. The reported pro- 
cedure exhibits good accuracy, within the typical uncertainty range of any experimental 
determination. The approach allows to extract all the information required for under- 
standing the noise performance of the device without any restriction on the statistics of 
the noise sources. The results show the relevant noise phenomena from a new angle. 

Keywords: Noise modelling, microwave transistor, stochastic integration, HEMT circuit 
model. 

1. Introduction 

The noise performance of an active device up to millimeter-wave frequencies plays 
a basic role in evaluating the suitability of the latter for low-noise applications in 
advanced telecommunication, nuclear instrumentation and radio astronomy applica- 
tions. This aspect is gaining increasing importance also due to the implementation 
of nanometer scale devices where a higher degree of interaction among the intrin- 
sic noise sources is expected to affect the device noise behavior pQ. So far, the 
noise analysis approaches and the relevant modelling tools employ the represen- 
tation based on a second order statistics for all the noise sources located inside a 
transistor |2). This allows to manage the analytical problems and leads to satisfac- 
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tory results within the 20 GHz range for device dimensions not aggressively scaled, 
as for example FET's with gate length not smaller than 0.2 /im. 

In this paper we propose a new approach to the problem starting from a general 
description of all noise sources and solving the stochastic integral equation associ- 
ated to the equivalent circuit of the microwave device ■ From the time waveform 
of the output voltage we obtain the noise figure F as a function of the operating 
frequency. The noise figure, which is the link between microscopic noise character- 
istics and macroscopic (i.e. measurable) parameters 0], is calculated by the ratio of 
the spectral noise power density produced at the output of the noisy device to the 
spectral noise power density produced at the output of the ideal noiseless device. 
More specifically, our calculation were performed for a noise source impedance of 
50 f2, but the procedure can be applied for arbitrary impedance value. We de- 
termine the device noise figure F50 at different discrete frequencies over a given 
operating range. Then, we compare the values of F50 calculated by means the pro- 
posed method with data referring to characterization of an HEMT transistor series 
(NE20283A, by NEC) from 6 to 18 GHz at a low noise bias point 5 . The proposed 
procedure exhibits fairly good accuracy, within the typical uncertainty range of any 
experimental determination. It is worthwhile to note that for Gaussian statistics 
the conventional CAD tools and our approach coincide. However for non-Gaussian 
noise sources, i.e. when it is important to know higher order statistics, our noise 
modelling is appropriate. In fact with our approach we can use different statistics 
of the noise sources and therefore obtain all the moments of the output voltage as a 
function of the cross-correlation functions of the noise sources. This technique can 
thus be applied to extract all the information required for a complete knowledge of 
the noise performance of the device without any restriction on the statistics of the 
noise sources, including also cross-correlated noise sources. 

2. HEMT Circuit Model 

High electron mobility transistors (HEMT) are GaAs or In-based field-effect tran- 
sistors for application at microwave and millimeter wave frequencies. These devices 
began to set records for noise performance since their introduction in the early '90s 
and therefore represent key elements in all low-noise circuits. The linear equiv- 
alent circuit which typically models the small-signal performance of HEMT's at 
microwave frequencies is shown in Fig. 1. This circuit represents the chip de- 
vice because access inductances and parasitic intererelectrode capacitances are not 
included. The values of the circuit elements are usually determined by computer- 
aided modelling procedures which may follow either optimization-based or direct 
extraction criteria |SJ[71|S]. More recently, approaches based on artificial neural 
networks have been successfully employed jHJ- We had previously characterized 
an HEMT transistor series (NE20283A, by NEC) over the 6-18 GHz frequency 
range at the low noise bias conditions suggested by the manufacturer (Vds = 2-0 V, 
Ids = 10 mA). By microwave measurements, we determined both the scattering 
and the noise parameters of the packaged device, i.e. the device sealed within a 
ceramic enclosure equipped with metallic leads. From the experimental data, we 
extracted a linear noisy model for the whole structure by a variable decomposition 
approach with a separate tuning for the channel noise source value modelled by 
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Fig 1. Small-signal circuit model of a chip HEMT with noise sources. 



an equivalent noise temperature assigned to the rd s resistor 9 . The values of the 
circuit elements and the noise temperature Td so obtained are listed in Tab. 1. 
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Table 1. Values of circuit elements of the chip HEMT model at Vds = 2.0 V and Ids = 10mA. 



All the other resistors in the model contribute with thermal noise typical of 
equilibrium conditions with device environment (i.e. S00K). The comparison be- 
tween the noise parameters measured and modelled for the packaged device data 
was presented in 0. We here report for the reader's convenience the plot of the 
noise figure in input matched condition in Fig. 2. To carry out our analysis we 
employ the model structure representing the chip device so eliminating the effect of 
the package found to be lossless (i.e. not containing noise sources). Among all noise 
parameters we use the F§q noise figure as information test to check the effectiveness 
of the procedure here presented. Therefore, the new results are compared with F$q 
data referring to the chip as simulated via CAD by use of the above circuit model. 



3. General Noise Analysis by Langevin Approach 

For the sake of simplicity and to validate our theoretical approach we suppose that 
all noise sources in the HEMT model are Gaussian white noise of thermal origin with 
zero mean, <5-correlated and variance per unit of frequency bandwidth equal to a 2 = 
AkT x R x , where k is the Boltzmann constant, R x the resistance and T x the equivalent 
temperature of the resistor. No correlation is assumed among the noise sources. The 
noise voltage generator e gen is the thermal noise associated with the real part of 
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Fig 2. Comparison between the measured noise figure F50 of the packaged HEMT transistor series 
NE20283A by NEC (o) and the F 50 simulated with CAD (-). 



the impedance Z gen = R gen + j u>X gen of the input termination at the reference 
temperature To = 290K, which is the standard temperature to determine the noise 
figure. The resistances r gs and r^s contribute to noise by assigning equivalent 
temperature T g and T^, respectively. The temperature of the intrinsic gate-source 
resistance r gs is assumed to be equal to the environment temperature (T g = To) [Hj- 
In order to determine the expression of the output noise voltage v Q as a function of 
the noise sources, the circuit network is analyzed in the Laplace domain. So that 
we have the Laplace transform of the zero-state (at t = 0) response of the output 
voltage 

6 

V (s) = Y,H i {s)E i {s), (1) 
i=i 

where s is a complex variable, Hi{s) are the transfer functions from the i-th noise 
generator to the output, and Ei(s) are the Laplace transforms of the voltage gen- 
erators with the following notation: 1 — > e gen , 2 — > e s , 3 — > e^, 4 — > e g , 5 — > e gs 
and 6 — > ed s - In the equivalent circuit of the HEMT, the small-signal g m transcon- 
ductance depends on the variable s as g m = g m oe~ ST , where g m o is the dc value 
of the transconductance and r is the transit time of the carriers to diffuse from 
source to drain. Since the value of r is a fraction of picosecond, and we are inter- 
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ested in investigating the frequency range below 20 GHz, where \s = iu>\ <C -, we 
approximate the exponential function as e~ ST ~ 1 — st + i(sr) 2 . Consequently the 
Hi(s) are rational functions of the complex variable s. By anti-transforming l|T]). 
we obtain the zero-state response of the output voltage in the time domain which 
assume the expression valid for t > 

6 pt 

v (t)=Y j hi(t-t')ei(t)df, (2) 

where the function hi(t) are the impulse responses corresponding to the i?i(s), 
which assume the expression 

hi{t) = h i0 6{t) + h ix e~~ + h i2 e~~ + h a eT~ , t > 0. (3) 

Here 5(t) is the Dirac <5-function. The time constants Tj of the impulse response 
and the coefficients hij depend on the circuit parameters and are obtained from the 
poles of the network functions if, (s) (see Appendix A). For the investigated model 
the following inequalities are satisfied 

t~3 <C tit 2 and uo <C — , (4) 

T3 

thus the relation © can be approximated as: 

hi(t) h lQ S(t) + h a e~~ + h l2 e~~ . (5) 

By using Eqs. (J2J and © the output voltage of the HEMT model in the time 
domain can be expressed as the following stochastic integral equation in the Ito 
sense 

6 pt 

v (t) = Y J [h i oe i (t)+ h. l {t-t')e l {t)dt l ] 1 (6) 
i=i ^ 

where the function h t contains only the last two terms of the expression ©.The 
Eq.(6) is the main result of our paper. It gives the stochastic output voltage of 
the HEMT device as a function of the different stochastic processes which represent 
the noise sources. By inserting the different statistics of the noise sources in Eq.(6) 
we can extract all the moments of the output voltage v (t) required for a complete 
knowledge of the noise performance when higher order statistics occurs in the output 
signal of the devices. As an example we report in Appendix B the expression of 
the second moment and the general expression of the nth moment of the output 
voltage as a function of the cross-correlation functions of the noise sources. As a 
first step of investigations and to validate our stochastic approach we consider all 
noise generators ei(t) present in the model as uncorrelated Gaussian white noises, 
with the usual statistical properties: {e%{t)) = and (ei(t) e,i{t )) = <rfS(t — t ). Here 
af is the variance per unit of frequency bandwidth of the i-th generator. Therefore 
the time waveform of the output voltage can be obtained by numerical integration 
of Eq. JHJ. In Eq. 10 the terms ei(t) dt = dWi(t) are Wiener processes associated 
with the device noise sources The numerical simulation of Eq.© is performed 
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by choosing the time step At = mm{ri, T2}/10, in order to prevent aliasing. At 
each temporal step t n — n At, by using the algorithm for Gaussian random variables 
[TT] . the Wiener increments are calculated as AWj(t„) = ai£(£„)vAt- The digital 
implementation of the impulse response JSJ allows to calculate the output voltage 
by the discrete-time convolution between the sequences h^tn) and the noise sources 

6 n 

gg^^-^AWift*). (7) 

i=l k=Q 

The duration of the sequence /ij(tn) increases whit n, therefore we consider only 
the first n max terms of the sequence h^tn) to save computation time. That is 

6 n 

Vo{t n ) = Y, [hioe l (t n ) + h' l {t n - k )AW l (t k )}. (8) 

The time discretization and the finite duration of the impulse response (FIR) in- 
troduce an error. In order to evaluate this error the z-transform of the sequence 
hi{to),hi{tx),hi{t2) hi(t nmax ) is considered 

"ma g 

H di (z)= hi(t n )z- n . (9) 

n=Q 

The error is evaluated by comparing the module of the Fourier transform of the 
continuous transfer function, i.e. 



ha 



Hi(iu) = h i0 + . + " a . , (10) 

1/Ti+lUJ 1/T2+10J 

with the Fourier transform of the sequence hi(t n ), i.e. H dl (z = e iujAt ), at different 
values of n max and for each generator e,-. The comparison between the continuous 
and the discrete transfer function is shown in Fig. 3, where we report the normalized 
Fourier transform of the transfer function versus the normalized frequency. In this 
figure we use n max = m x max{ri,T2}, with m integer, and f s — 1/(2 At) is 
the Nyquist frequency. As m increases, the difference between continuous and the 
approximate impulse response becomes negligible. At m = 6 the error introduced 
by the FIR impulse response are less than 1%. We use this value of m in our 
simulations. 



4. Results 

The samples of the output noise voltage waveforms were calculated for an input 
termination 50 from Eq. © by using C-language routines. To validate the 
approach proposed in this paper we compare the theoretical results of the noise 
figure F50 with that obtained by CAD analysis in the 6-18 GHz frequency range. 
The temporal noise sequence was generated by using a time step At of 5.57 x 10~ 13 
sec and Gaussian random number generator. The Nyquist frequency is 0.898 THz, 
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Fig 3. Normalized module of the Fourier transform of the transfer function from the input generator 
to voltage output for the continuous and discrete impulse response and for different values of the 
integer m. In the x-axis the frequency is normalized with respect to the Nyquist frequency and 
the values are restricted to the 6 — 18 GHz frequency band. 



thus N = 10 5 points are needed at least to investigate the microwave frequency 
range. The noise figure is calculated by using the following relation |12|: 

r _ Pon(f) , s 

^50 777 yi-i-i 

where p n(f) is the power spectral density of the output voltage at the frequency / 
when all the noise sources are activated and p f / (/) is the power spectral density 
obtained from the noise source of the input termination only, both in V 2 /Hz units. 
In order to improve the estimation of the spectral power density given by Ijllfl at 
a given value of frequency /, we use the mean value of the noise power evaluated 
within the frequency bandwidth BW around / 



rf+BW/2 

<?*(/)> = -5TJ7 / PM)df- (12) 

lf-BW/2 



BW 



The output noise waveforms of the HEMT model, with (v on ) and without (v ff) 
the internal noise sources, over a sample time interval of 20 ns are reported in Fig. 
4. The total time period used for calculation is N At, where N is usually chosen to 
be 1 x 10 5 . 

By using BW = 90 MHz for noise power estimation, we calculate the F50 as a 
function of the operating frequency in the 6-18 GHz range. In Fig. 5 we compare the 
noise figure behavior obtained with our approach with that obtained by computer- 
aided analysis of the model. 
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Fig 5. Noise figure F50 calculated with Langevin approach and by using BW = 90 MHz for noise 
power estimation ( — ) compared with F50 obtained by computer-aided analysis of the model ( — ) . 
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Good agreement is exhibited by the two noise figure curves. We use 100 real- 
izations of time waveforms to determine each noise figure behavior over the given 
frequency range. The results obtained with the stochastic approach may be im- 
proved by averaging over a greater number of realizations but with more computer 
time consuming. We note that the noise figure behavior obtained with the stochas- 
tic approach could be thought as the output of a noise figure measuring instrument 
where averaging and smoothing options are available to minimize display jitter. The 
procedure here described can be used with any input termination value to obtain 
more than one noise figure value. By choosing four different termination values, 
the extraction of all the standard noise parameters can be performed by using the 
equation 

F ^' F ° +ir " nJFii-iLp) (13) 

where T gen is the generic termination value expressed in a reflection coefficient form 
and Fq, T (magnitude and angle) and r n are the four noise parameters usually 
determined by complex experimental procedures |lri| . 

By calculating the noise parameters vs. frequency, a complete knowledge of the 
noise behavior of any two-port is obtained from a circuit point of view. The noise 
parameters can be employed either for extraction of device intrinsic noise sources 
or for CAD of low-noise circuits. 

5. Conclusions 

We have presented a noise analysis procedure for microwave transistors. This pro- 
cedure is based on Langevin approach and takes the general statistical properties of 
the internal noise sources into account. The device is represented by its equivalent 
circuit with noise generators to a chip level. The noise figure at a given input ter- 
mination is determined by first integrating the noise voltage equations in the time 
domain and by transforming them into spectral power densities over the frequency 
range of interest. The results are compared with the simulated data, obtained with 
a commercial CAD, of the chip model of an HEMT device series (NE20283A, by 
NEC) from 6 to 18 GHz at a low noise bias point. The proposed procedure exhibits 
good accuracy within the typical uncertainty range of any experimental determi- 
nation. Our Langevin approach gives the same results of noise analysis approaches 
based on second order statistics. But when higher order statistics becomes rele- 
vant to the noise analysis our approach is suitable. Moreover the cross-correlation 
functions between the noise sources of the device are explicitly included in our the- 
oretical approach (see Appendix B), while in conventional CAD methods their are 
inserted in a heuristic way. Such a general approach can be extended to analyze 
noise in other devices without restrictions as to the frequency range, by specify- 
ing the device equivalent circuit and the statistical properties of the internal noise 
sources. 
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Appendix A. 

Because of the topological complexity of the HEMT model network it is not possible 
to obtain the expression of the time constants Tj in a closed form. Indeed the Tj are 
obtained from the poles of the transfer function Hi(s), which are the roots of the 
third-order algebraic equation 

2(1 + (Cg s r gs + Cgdfds + CdsTds + Cg S Rg + CgdRg + RgCTlCgd + RgenCgs 
+gm,OrdsCgdRg + R S Cg S + QmoRsfdsCgd + R S Cgd + 



Principato et al. 



+ (C 'dsT dsC 'g s r g S + C dsT dsC gdRg + CgdTdsCgsRg + Cds r dsCg S Rg + C g S T g S C gdRg 
~9mO r ds T CgdRg + T dsC 'g 3 r g S C ' gd + R ds^ dsC gd + Rs^dsCg S Cgd + RsCds r dsCg S 

^rTgsRsCgsCgd — gmoRs r ds T Cgd + CgdTdsRgenCgs + C ds r ds RgenC gs 
+ C <j«F dsRgenC 'gd + C g S T gs C gdRgen ~ 9mO r ds T C gdRg en) S 2 ) 
+ (Cgdrds ( 9mO T + ^TgsGgsCds) [Rg + + Rgen)) S — 0. 

(A.l) 

Appendix B. 

In the HEMT model the only plausible cross-correlated noise sources are the gate- 
source e ffs and the drain-source e^s [H!< Do to this hypothesis we obtain from Eq. 
(6) the general expression for the second moment of the output voltage Vo(t) as a 
function of the cross-correlation function between the e gs and the e^s noise sources 
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(B.2) 

where we use the property of the Dirac (^-function. The general expression of the 
moments of Vo(t) with the same above-mentioned hypothesis is 
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